AO - A 1 29  395 
UNCLASSIFIED 


FEASIBILITY  OF  APPLYING  THE  FINITE  ELEMENT  ADAPTIVE  1/  | 

RESEARCH  SOLVER  (FEAR. . (U)  MARYLAND  UNIV  COLLEGE  PARK 
INST  FOR  PHYSICAL  SCIENCE  AND  TECH..  0  A  GIGNAC  ET  AL . 

MAY  83  DTNSRDC/CMLD-83/11  NOO 167-82-M-0743  F/G  9/2 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-  1963-A 


FEASIBILITY  OF  APPLYING  THE  FEARS  PROGRAM  TO  THE  PLATE  BENDING  PROBLEM  DTNSRDC/CMLD-83/11 


( - 

DAVID  W.  TAYLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CENTER 

Bath* ids,  Maryland  20084 


m 

co 

<N 


! 

FEASIBILITY  OF  APPLYING  THE  FEARS  PROGRAM 
TO  THE  PLATE  BENDING  PROBLEM 


by 


Donald  A.  Gignac 
David  W.  Taylor  Naval  Ship  R&D 


and 

Ivo  Babuska 

University  of  Maryland 
Contract  N00167-82-M-0743 


Center 


OTIC 

ELECTE' 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


COMPUTATION,  MATHEMATICS,  AND  LOGISTICS  DEPARTMENT 
DEPARTMENTAL  REPORT 


May  1983 


DTNSRDC/CMLD-83/  11 


J 


NOW  OTNSROC  5*0230  (2-80) 
liupwiwM  3900/46) 


83  06  16  05  5 


MAJOR  DTNSRDC  ORGANIZATIONAL  COMPONENTS 


GPO  8*»  H? 


NOW-OTNSRDC  MO 2/21  (2-MI 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FACE  (V»«i  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


r.i  i.M.itf’Tiri  1 1. 


DTNSRDC/CMLD-83/11 


TITLE  (end  Subtitle) 

FEASIBILITY  OF  APPLYING  THE  FEARS  PROGRAM 
TO  THE  PLATE  BENDING  PROBLEM 


.  AUTMONfaj 

Donald  A.  Glgnac 
I vo  Babuska 


s.  rcrformino  oroaniiation  name  ano  address 
David  Taylor  Naval  Ship  Research  and 
Development  Center 
Bethesda,  Maryland  20084 


n.  controlling  office  name  ano  address 


read  nmncnoMi 
BEFORE  COMPLETING  FORM 


REORIENT  S  CATALOO  NUMBER 


'•WHW' 


S.  TYRE  OF  RERORT  •  RERIOO  COVERED 


S.  RERFORMINO  ORO.  RERORT  NUMBER 


I 


MONITORING  AGENCY  NAME  A  ADDRESSflf  dltleront  from  Controlling  Olllco) 


RROORAM  ELEMENT.  PROJECT.  TASK 

Pr  o|  r  am  emen  E  2 
Task  Area  ZR0140201 
Work  Unit  1844-140 


12.  RERORT  DATE 

May  1983 


IS.  NUMBER  OF  PAGES 

27 


IS.  SECURITY  CLASS,  (ol  thlo  report) 

UNCLASSIFIED 


is«.  declassification/downgrading 
SCHEDULE 


l«.  DISTRIBUTION  STATEMENT  fol  thlo  Report) 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


17.  DISTRIBUTION  STATEMENT  (ol  Uio  mbmtrmct  entered  In  Block  30,  II  dlllerent  frooi  Roport) 


IS.  KCY  WORDS  (Continue  on  rereroo  olde  If  nocoooory  end  Identity  by  block  number) 


Finite  Element  Method 
Plate  Bending  Problem 
FEARS  Program 
Mesh  Refinement 
Error  Estimation 


Adaptive  Meshing 
Lame'  Equations 
Splitting  Method 
Biharmonlc  Equation 
Moments 


Shear  Forces 


20.  ABSTRACT  (Continue  on  teoeroe  olde  II  nocoooory  end  Identity  by  block  number) 

-^This  report  documents  an  Investigation  of  the  feasibility  of  the 
application  of  the  Finite  Element  Adaptive  Research  Solver  (FEARS)  computer 
program  to  the  Plate  Bending  Problem.  Two  methods  of  reducing  this  bihar¬ 
monic  problem  to  an  elliptic  system  of  two  second  order  partial  differential 
equations  are  considered.  The  first  is  the  splitting  method  and  the  second 
is  the  transformation  to  the  Lame'  system  of  elasticity  equations.  The  FEARS 
program  is  used  to  solve  these  reduced  systems  for  three  examples.  The 

(continued  on  reverse  side) 


FORM 
I  J  AN  71 


COITION  OF  I  NOV  ••  It  OSIOLtTB 
S/N  01 02-  L  F  -01 4-6601 


_ UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIt  RAOE  (WA**  I 


r 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  RAOC  f*hm  Data  Balara<0 


(Block  20  continued) 


\ 


esults  are  analyzed  and  discussed  with  regard  to  the . computation  of 
displacements,  moment,  and  shear  forces.  While  the  Lame*  system 
approach  is  well-founded  theoretically,  it  poses  some  computational 
problems  with  regard  to  accuracy.  The  authors  feel  that  the  splitting 
method  when  monitored  by  error  indicators  is  the  preferred  method  for 
FEARS  users  even  though  the  theory  of  this  method  has  not  been  com¬ 
pletely  developed  yet. 


* 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF“THIS  FASCISM*  ®a*a 


TABLE  OF  CONTENTS 


Page 


LIST  OF  FIGURES . iii 

LIST  OF  TABLES .  iii 

ABSTRACT  .  1 

ADMINISTRATIVE  INFORMATION  .  1 

INTRODUCTION  .  1 

ADAPTIVE  FINITE  ELEMENT  SOLVER  FEARS  .  2 

SAMPLE  PROBLEM  .  4 

REDUCTION  TO  THE  LAME*  SYSTEM  OF  ELASTICITY  EQUATIONS  .  5 

SPLITTING  METHOD . ' .  8 

NUMERICAL  EXPERIMENTATION  .  11 

REDUCTION  TO  THE  LAME'  SYSTEM  OF  ELASTICITY  EQUATIONS  .  13 

SPLITTING  METHOD  .  14 

CONCLUSIONS .  19 

ACKNOWLEDGMENTS  .  19 

REFERENCES . 21 

LIST  OF  FIGURES 

1  -  Partition  of  Domain  ft  1 .  11 

2  -  Domain  ....  . . 12 

3  -  Partition  of  Domain  .  12 

LIST  OF  TABLES 


1  -  The  Value  of  the  Second  Partial  Derivatives  of  Z  and  W  at  the 


Point  (.53125,  .53125)  (Example  3,  Z  =  Z(1)) .  14 

2  -  Example  1:  Uniform  Mesh . 15 

3  -  Example  2:  Uniform  Mesh .  16 

4  -  Example  3:  Uniform  Mesh .  17 

5  -  Example  3:  The  Adaptive  Mesh  . .  18 


iii 


ABSTRACT 


This  report  documents  an  investigation  of 
the  feasibility  of  the  application  of  the  Finite 
Element  Adaptive  Research  Solver  (FEARS)  computer 
program  to  the  Plate  Bending  Problem.  Two  methods 
of  reducing  this  biharmonic  problem  to  an  elliptic 
system  of  two  second  order  partial  differential 
equations  are  considered.  The  first  is  the  split-  j 
ting  method  and  the  second  is  the  transformation  to 
the  Lame*  system  of  elasticity  equations.  The  FEARS 
program  is  used  to  solve  these  reduced  systems  for 
three  examples.  The  results  are  analyzed  and  dis¬ 
cussed  with  regard  to  the  computation  of  displace¬ 
ments,  moment,  and  shear  forces.  While  the  Lame' 
system  approach  is  well-founded  theoretically,  it 
poses  some  computational  problems  with  regard  to 
accuracy.  The  authors  feel  that  the  splitting 
method  when  monitored  by  error  indicators  is  the 
preferred  method  for  FEARS  users  even  though  the  r. 
theory  of  this  method  has  not  been  completely 
developed  yet.  j 
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This  work  was  performed  under  David  W.  Taylor  Naval  Ship  R&D  Center's 
(DTNSRDC)  Independent  Research  IVogram,  Program  Element  61152,  Task  Area 
ZR014020 1,  DTNSRDC  Work  Unit  1844-140.  Professor  Ivo  Babuska  of  the  Institute 


for  Riysical  Science  and  Technology,  University  of  Maryland  (Contract 
N00167-82-M-0743)  ia  a  contractor  to  DTNSRDC. 


INTRODUCTION 

The  Finite  Element  Adaptive  Research  Solver  (FEARS)  computer  program  i3 

1  2* 

an  adaptive  finite  element  solver  with  a  posteriori  error  estimates.  *  It 
solves  boundary  value  problems  for  an  elliptic  system  of  two  partial  differ¬ 
ential  equations  with  two  independent  variables,  using  first  degree  elements 
which  have  continuous  displacements  across  element  boundaries.  Elements  with 
such  continuous  displacements  are  said  to  be  (?. 

The  classical  plate  bending  problem  is  described  by  a  partial  differential 
equation  of  the  fourth  order ,  the  biharmonic  equation.  It  should  be  possible 
to  transform  this  biharmonic  problem  into  an  elliptic  system  of  two  second 

*  A  complete  list  of  references  is  given  on  page  21. 
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order  partial  differential  equations,  but  the  problem  of  obtaining  a  family  of 

such  transformations  has  not  been  solved  yet  for  the  general  case. 

The  biharmonic  problem  can  be  solved  through  a  reduction  to  a  system  of 

second  order  equations  by  various  methods:  the  Reissner-Mindlin  refined  plate 
3  4 

formulation  which  provides  three  fields  of  variables  (displacements  and 

two  rotations),  the  Hermann-Mioshi  method3’  ’’*  ’ ’  and  the  Hermann-Johnson 

method'’’*’’ 10  which  provides  a  system  of  four  equations.  Other  methods  of 

solving  biharmonic  equations  are  described  by  Scholz. 11  In  addition,  the 

12  13 

splitting  method  first  considered  by  Glowinski  and  Mercier  and  further 
8  9  11  14  15 

developed  by  others  ’  '  provides  a  system  of  two  second  order  equa¬ 

tions.  Finally  the  biharmonic  problem  can  be  transformed  into  the  Lame*  system 
of  elasticity  equations  by  exploiting  the  fact  that  the  Airy  function  is  bihar¬ 
monic.  The  last  two  approaches  provide  a  boundary  value  problem  solvable  by 
the  FEARS  program. 

This  report  studies  the  performance  of  these  two  last  approaches  with 
respect  to  the  FEARS  program  and  assesses  their  advantages  and  disadvantages. 
The  second  section  contains  basic  information  about  the  FEARS  program.  The 
sample  problem,  the  method  of  transformation  to  the  Lame'  equations,  and  the 
splitting  method  are  described  in  the  third  section.  The  numerical  experi¬ 
ments  are  described  and  some  conclusions  and  recommendations  are  presented  in 
the  fourth  and  fifth  sections  respectively. 

ADAPTIVE  FINITE  ELEMENT  SOLVER  FEARS 

1  2 

The  FEARS  program  has  been  described  in  two  informal  reports.  ’  It 
solves  boundary  value  problems  for  linear  self-adjoint  elliptic  systems  of 
partial  differential  equations  of  second  order  with  two  field  variables  in  two 
dimensions.  It  also  provides  error  estimates  for  the  finite  element  solution 
with  respect  to  the  unknown  exact  solution  of  the  problem. 

FEARS  deals  with  domains  which  are  the  union  of  the  interiors  of  a  small 
number  of  open  curvilinear  rectangles  bounded  by  circular  arcs  or  straight  line 
segments.  These  curvilinear  rectangles  are  called  2-D  domains  and  are  denoted 
by  D^,  J=1,...,N2.  The  open  arcs  are  called  1-D  domains  and  are  denoted  by  D* , 
j*1,...,Nj.  The  vertices  of  the  2-D  domains  are  called  0-D  domains  and  are 
denoted  by  Dj,  j=i . Nq. 
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The  FEARS  prograa  solves  problems  described  in  the  weak  form 


l  J  +  i  “V-Z 


+  £  {  e.(s)Vds 

i  J 1  1 

D. 


where 


and 


3U  - 
9Z 


V  an<*  Uz  are  defined  similarly. 


(1) 
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Further 


^  is  a  4  t  4  symmetric  (constant)  matrix, 

Bi  is  a  2  x  4  (constant)  matrix, 

is  a  2  x  2  symmetric  (constant)  matrix, 

Yj  is  a  2  x  2  symmetric  (constant)  matrix, 

Di^x1,x2^  ^  a  **  x  1  vector  valued  function  with  components  S^x^Xg) 

(k  =  1,2, 3,4), 

Ei*x1'x2*  is  8  2  X  1  vector  valued  function, 

e1(s)  (where  s  is  arc  length)  is  a  1  x  2  vector  valued  function. 

The  trial  space  used  for  U  and  the  test  space  used  for  V  allow  the  imposi¬ 
tion  of  both  essential  and  natural  boundary  conditions.  For  a  discussion  of 

the  physical  significance  of  the  matrices  with  examples,  the  reader  is  referred 

2 

to  the  User's  Guide. 

The  admissible  error  norm  for  the  a  posteriori  estimate  of  the  error 

•  s  u_.-u  (where  u__  and  u  denote  the  finite  element  and  exact  solutions, 

FE  o  FE  o 

respectively)  is  j | | e| | j  ^  where 

IIWII*-  l2> 

Di 

FEARS  uses  C°  first  degree  elements  of  bilinear  type  on  curvilinear  rec¬ 
tangles  (which  are  mapped  into  bilinear  elements  on  the  master  square)  and 
adaptively  constructs  the  meshes  by  refinement  where  appropriate. 

SAMPLE  PROBLEM 

Let  OCR2  be  a  bounded  simply  connected  region  with  a  boundary  3£1. 

The  problem  is  to  find  a  function  U  satisfying  the  biharmonic  equation 

AAU  =  f  on  n  (3) 


and  the  boundary  conditions 


on  30 


(4) 
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Assum  that  there  ia  a  saooth  function  Z  such  that  g=Z  and  h  =  —  on  an. 

dn 

U  is  the  deflection  of  the  plate.  The  second  derivatives  of  U  are  the  bending 
aaaents  which  are  of  special  interest.  At  present  only  the  problem  of  a 
built-in  plate  with  boundary  conditions  of  the  type  given  by  Equation  (4)  will 
be  considered. 

REDUCTION  TO  THE  LAME'  SYSTEM  OF  ELASTICITY  EQUATIONS 

Assiae  that  fsO  in  Equation  (3).  Then  U  is  a  biharoonic  function  which 
can  be  regarded  as  the  Airy  function  for  the  plane  elasticity  problem.  There¬ 
fore,  there  exist  functions  u^Ug  defined  on  0  such  that 


and  the  functions  u1§u2  also  satisfy  the  well-known  Lame’  equations  (for  v  =  o 
where  v  is  Poisson's  ratio). 

The  functions  g  and  h  must  be  transformed  into  the  boundary  conditions  for 
Uj.Ug.  It  is  not  hard  to  show  that  u^Ug  satisfy  the  traction  condition  on 
the  boundary.  However,  it  must  be  emphasized  that  Q  is  assumed  to  be  simply 
connected.  (If  the  plate  were  free,  then  the  transformed  problem  would  provide 
the  prescribed  displacements  for  u^.Ug  on  3D.) 

The  relationship  between  the  traction  and  the  functions  g  and  h  is  not 
complicated,  but  it  involves  differentiation  along  the  boundary.  The  assump¬ 
tion  that  the  functions  g  and  h  are  defined  in  terms  of  the  function  Z  simpli¬ 
fies  the  transformation .  The  framework  of  the  FEARS  program  is  especially  well 
suited  to  take  advantage  of  this  assumption. 

The  weak  formulation  of  the  problem  for  the  functions  u^Ug  is  a  standard 
one.  The  functions  Uj.u^  satisfy  the  variational  form  of  the  problem  (virtual 
work  condition)  and  are  not  subject  to  any  essential  boundary  conditions. 


More  precisely,  u..,u_e  H  (where  H1  is  the  usual  Sobolev  space)  are  functions 

1  £  i 

such  that,  for  Vj,v2  eH  , 


B1(U1’U2:  V1,V2>  =  F(vi*v2)  (6> 

Here  the  bilinear  form  B  is  the  usual  form  of  virtual  work  obtained  by  taking 
E=1,  v  —0 •  (E  is  Young’s  modulus  of  elasticity  and  u  is  ftoisson’s  ratio.)  Then 


VU1*°2 


+ 


;  vrv2) 


3u,\/<*v 


■)(; 


dX)dx2 


(7) 


The  functional  FCv^Vg)  expresses  the  virtual  work  of  the  traction  forces 
on  the  boundary.  Integration  by  parts  provides  FCv.j.v^  in  terms  of  the 
function  Z.  Thus 


FCv^Vj) 


ff  3V1  32Z  9v2  32Z 

V  K  3x2  "  3xl  3X1X 


3v 

3x 


1  32Z 


2  Sx^x 


2 


+ 


3V2  32z1 


dx^dx2 


(8) 


This  functional  is  of  a  form  which  the  FEARS  program  can  handle  if 


.2 
3  Z 


Bx^x 


2 
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The  conditions  given  by  Equations  (6),  (7).  and  (8)  determine  u.|,u2 
uniquely  only  up  to  rigid  body  motion.  Accordingly,  a  displacement  is  speci¬ 
fied  at  one  point  and  a  rotation  through  another  point. 

When  the  functions  u^u^  have  been  obtained,  the  second  derivatives  of  U 
(the  bending  moments)  are  provided  by  Equation  (5).  An  additional  integra¬ 
tion  must  be  performed  to  obtain  the  values  of  U.  If  however,  the  bending 
moments  are  the  main  concern,  then  the  functions  u^,u2  anc'  Equation  (5) 
provide  them  directly. 

So  far  it  has  been  assumed  fsO.  Now  assume  that  f$0  and  g=h=0.  In  this 
case  a  particular  solution  is  needed,  which  often  (e.g.,  for  f=1)  can  be  easily 
found.  If  this  particular  solution  is  denoted  by  Z,  then 

U  =  Z  -  W 


where 


AAW  =  0  on  ft 


and 


1  =  W*  TT  "  7T  on  3£2* 
3n  3n 


This  is  exactly  the  case  discussed  earlier.  Now 


2 

3*U 


32Z 


32w 

Sxjaxj  ’ 


i,j  ■  0,1,2  such  that  i+j  *  2 


(9) 


This  approach  has  the  following  advantages  and  disadvantages: 

1)  Advantages.  The  finite  element  method  for  solving  elasticity  pro¬ 
blems  is  well-founded  theoretically.  The  error  with  respect  to  the  energy 
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norm  (which  is  essentially  the  I*2  norm  for  the  bending  moments)  can  be  reli¬ 
ably  estimated.  The  postprocessing  approach  allows  the  computation  of  the 
moments  at  selected  points  with  generally  better  accuracy  along  with  an  error 
estimate.  Likewise  the  stress  intensity  factors  can  also  be  computed. 

2)  Disadvantages.  If  f$0,  then  a  particular  solution  Z  must  be  known. 
Obtaining  the  particular  solution  Z  for  a  general  f  could  be  very  laborious. 
However  obtaining  Z  for  a  particular  f  in  most  applied  problems  is  not  that 
difficult.  Often  Z  can  be  found  by  inspection  or  some  artifice.  Similarly 
the  computation  of  the  displacements  (the  values  of  U)  could  also  be  quite 
laborious.  In  addition,  the  accuracy  of  the  numerical  solution  depends  on 
the  selection  of  the  particular  solution.  Since  values  of  the  moments  are 
obtained  from  Equation  (9)  by  subtracting  numbers  of  approximately  the  same 
magnitude,  the  relative  error  in  the  computation  of  W  can  result  in  a  much 


larger  relative  error 


Accordingly,  it  is  necessary  to  compute 


W  with  greater  relative  accuracy.  If  the  domain  flis  not  simply  connected, 


then  additional  difficulties  arise. 


SPLITTING  METHOD 

Although  the  theory  of  this  method  has  been  studied  by  various 
authors,  ^  1  the  method  is  still  not  very  well  understood.  Today's 

theoretical  results  presuppose  assumptions  which  are  rather  restrictive  and  it 
is  not  clear  whether  these  assumptions  are  really  necessary.  These  basic 
assumptions  are 


( 1 )  The  elements  used  are  at  least  of  degree  2. 

(2)  The  mesh  is  quasi-uniform. 

(3)  The  domain  0  is  convex. 

However,  these  assumptions  do  not  seem  necessary.  Computational  results 
obtained  in  the  course  of  this  investigation  strongly  suggest  that  the  assump¬ 
tions  are  sufficient  and  not  necessary  conditions. 


For  simplicity,  consider  the  problea  given  by  Equations  (3),  (4)  with 
gshsO.  The  method  consists  in  writing  Equation  (3)  in  the  fora 


Au 


1 


(10) 


where  the  boundary  conditions  are 


u 


1 


s  0 


(11) 


(12) 


The  bilinear  fora  is 


B2*u1,u2;  vl *v2^  * 

ff/3ul  3V2  3ul  3v2 
Ji\3x1  SXj  +  3x2  3x2  +  u2v2 


3u 

lx 


,  3v.  3u,  3v  \ 

—  — -  +  — -  — L  I  dx  dx„ 
1  3xi  3x2  3x2/  12 


(13) 


and  the  functional  is 


F2<  W  *  -//fv i  dxidx; 
n 


(14) 


Then  the  problea  can  be  written  in  the  standard  weak  fora 

B2^U1,U2J  V1,V2^  "  F2^V1,V2^ 


(15) 


9 


with 


u^,v1  e  H1,  Ug.Vg  e  H1  where  H1  =  |u  e  H1  |  u  =  0  on  3i)f 

The  boundary  condition  given  by  Equation  (11)  is  essential  and  the  boundary 
condition  given  by  Equation  (12)  is  natural.  The  bilinear  form  B.  is  not 
positive  definite  and  Equation  (15)  is  a  "saddle  point"  problem. 

The  bilinear  form  and  the  functional  F 2  can  be  used  by  the  FEARS  pro¬ 
gram.  However,  since  FEARS  employs  elements  of  degree  1,  the  first  condition 
is  violated.  The  numerical  examples  of  this  investigation  also  violate  the 
second  and  third  conditions. 

There  are  various  options  for  error  norms  in  the  FEARS  program.  The 

norms 


can  be  used.  For  a  nonconvex  domain  the  norm  |  [ u | | ^  of  the  exact  solution  u 
is  infinite.  If  the  solution  is  smooth  then  either  of  the  norms  ||.||  or 
||.||2  «ay  be  used.  If  the  solution  has  singular  behavior,  then  only  the 
norm  |  |.|  ]2  can  be  used. 

This  approach  has  the  following  several  advantages  and  disadvantages: 

1)  Advantages.  The  problem  can  be  solved  by  FEARS  without  obtaining 
a  particular  solution.  This  approach  avoids  the  difficulty  stemming  from  the 
arbitrary  selection  of  a  particular  solution  and  it  also  avoids  the  problem  of 
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significance  due  to  subtraction  encountered  with  the  previous  nethod .  The 
aethod  directly  coaputes  the  displaceaent  and  its  Laplacian.  The  assiaption 
that  the  doaain  is  siaply  connected  is  not  required  for  this  aethod. 

2)  Disadvantages.  The  aethod  is  not  coapletely  understood  theoreti¬ 
cally.  The  saddle  point  approach  does  not  guarantee  aonotonic  convergence 
(which  occurs  when  the  bilinear  fora  is  positive  definite).  At  present  the 
convergence  of  the  aethod  has  not  been  established  under  the  conditions  of  the 
FEARS  coaputation.  The  second  derivatives  (which  are  the  aoaents)  oust  be 
computed  by  postprocessing.  The  error  norms  used  by  FEARS  do  not  monitor  the 
accuracy  of  the  moaents. 


NUMERICAL  EXPERIMENTATION 

Consider  the  three  examples  given  by  Equations  (3)  and  (4)  where  f=64  on 
domains  and  g=h=0  on  their  boundaries  3^,  d(l2,  a^. 

Example  1  -  Uniformly  Loaded  Circular  Plate 

«1  =  j  (xrx2)  |  xt  ♦  x2  <  l} 

ft  1  is  partitioned  into  the  five  2-D  domains  shown  in  Figure  1. 


Figure  1  -  Partition  of  Domain  0 1 
The  exact  solution  is  known  to  be 

U  »  (1-(x1  «.  x2))2 


11 


Example  2  -  Uniformly  Loaded  Square  Plate 


fl2  =  {(x1fx2)  |  |  xA  |  <  1.  1  *  1,2} 

Figure  2  shows  the  unpartitioned  domain  ftg. 


The  solution  is  not  available  in  closed  form  though  there  are  analytical 
pressions  for  U  in  series  form. 


Example  3  -  Uniformly  Loaded  Square  Plate  With  a  Slit 

n3  =  a2  -  { (xrx2)  |  -1  <  x1  <  0,  x2  =  0} 

3  is  partitioned  into  four  equal  squares  as  shown  in  Figure  3* 


The  exact  solution  is  not  available  for  Example  3. 


The  exact  solution  is  very  smooth  for  Example  1,  fairly  smooth  at  the  cor¬ 
ners  for  Example  2,  and  has  a  strong  singularity  at  the  origin  for  Example  3. 
The  trial  and  test  functions  are  bilinear  for  Examples  2  and  3t  and  for 
Example  1  they  are  mapped  onto  bilinear  functions.  These  examples  have 
different  axes  of  symmetry  which  were  utilized  to  a  limited  extent,  but  the 
material  in  the  tables  is  organized  as  though  no  symmetries  were  present. 

The  two  approaches  of  the  third  section  will  now  be  discussed. 


REDUCTION  TO  THE  LAME'  EQUATIONS  OF  ELASTICITY 
We  use  the  following  functions  Z  and  S^ . 


‘"’■ft*  <Y 


sj15  =  4x2  ♦  12x2 


S(1)  -  S(1)  -  -A*  . 

S2  "  S3  -  “®*1*2 

S^}  =  12x2  ♦  4x2 


Z(2)  =  8/3xJ* 


S<2)  =  s2(2)  =  S<2> 


=  0 


.(2) 


=  32 


Z(3)  =  8x2x2 


,(3> 


16  x‘ 


.(3) 


-32x^2 


, (3)  , 


=  16X2 
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The  Method  provides  the  expected  results.  The  selection  of  Z  barely  influenced 
the  relative  accuracy  of  the  computed  function  W  with  respect  to  the  energy 
norm  when  the  same  number  of  elements  was  used.  The  error  estimation  was  very 
effective.  The  main  problem  of  the  method  is  related  to  the  subtraction  of  the 
values  of  Z  and  W.  Table  1  shows  typical  results  for  Example  3  (Z  =  Z^) 
with  a  mesh  of  1024  elements.  This  table  displays  the  values  of  the  second 
derivatives  of  the  functions  Z  and  W  at  the  point  (.53125.  .53125). 

TABLE  1  -  THE  VALUES  OF  THE  SECOND  PARTIAL  DERIVATIVES  OF  Z  AND  W 
AT  THE  POINT  (.53125,  .53125)  (EXAMPLE  3,  Z  =  Z(1>) 

_3i  j£  a2 

2  2  — - - 

3XJ  3x2  axj  3x2 

W  5.89848  5.72560  1.81340 

Z  4.51562  4.51562  2.25781 

31 2 3U 

Now  an  accuracy  of  the  order  of  5%  in  -  ,  where  U=Z-W,  means  that  W  must 

i  j 

have  an  accuracy  in  the  neighborhood  of  1-21.  This  restriction  on  W  could  be 
prohibitive.  This  problem  arises  in  the  most  important  case  when  f|o.  If  f=0, 
then  this  problem  is  usually  not  so  severe. 

SPLITTING  METHOD 

As  noted  in  the  previous  section,  convergence  is  not  theoretically  guaran¬ 
teed  for  this  method.  The  solution  was  computed  using  both  a  uniform  mesh  and 
an  adaptive  mesh  with  respect  to  the  error  estimator  of  the  |  |  •  |  |  2  nonD  given 
by  Equation  (17).  The  results  are  displayed  in  Tables  2-5.  In  these  tables 
the  computed  values  of  u^  and  u^  on  the  line  0A  are  given.  (As  the  meshes  are 
refined,  more  values  on  0A  are  computed.)  The  headings  of  these  tables  are  as 
follows : 

(1)  The  x.|-  coordinate  of  the  point  (x^,0) 

(2)  The  computed  value  of  u^x^.O) 

(3)  The  computed  value  of  “2<*i,0) 
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TABLE  2  -  EXAMPLE  1 :  UNIFORM  MESH 


A 


TABLE  3  -  EXAMPLE  2 :  UNIFORM  MESH 


N 

_ $4 _ 1 

256 

1024 

X1 

Ujtx^O) 

u2(xltO) 

u1(x1»0) 

u2<xi°> 

Uj^XjO) 

u1(x1,0) 

0 

1.35644 

-9.61960 

1.31121 

-9.17025 

1.29959 

-9.05834 

.0625 

1.29076 

-8.99565 

.1250 

1.27571 

-8.91722 

1.26440 

-8.80644 

.1875 

1.22086 

-8.48732 

.2500 

1.21137 

-8.57152 

1.17114 

-8.13994 

1.16078 

-8.03261 

.3125 

1.06506 

-7.43436 

.3750 

1.00384 

-6.78387 

.994969 

-6.68236 

.4375 

.892184 

-5.76413 

.5000 

.812358 

-5.13491 

.785777 

-4.75807 

.778854 

-4.66495 

.5625 

.657709 

-3.36783 

.6250 

.536873 

-1.93536 

.532167 

-1.85355 

.6875 

.406453 

-1.00715 

.7500 

.297489 

1.57309 

.288256 

1.84719 

.285755 

1.91421 

.8125 

.176378 

4.21680 

.8750 

.0866609 

6.78649 

.0859307 

6.83454 

.9375 

.0235262 

9.79657 

1,0000 

0 

13.0037 

0 

13.1088 

0 

13.1334 
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TABLE  4  -  EXAMPLE  3:  UNIFORM  MESH 


N 

64  "1. 

256  j 

I  1024 

*1 

ul*xl°* 

u2(*!,0) 

UjUjO) 

u2(x1,0) 

u9(xj ,0) 

0 

0 

14.0294 

0 

20.5651 

0 

29.0652 

.0625 

.0366752 

9.02679 

.1250 

.097985 

3.  21982 

.0856739 

4.39889 

.1875 

.133943 

1.45009 

.2500 

.199010 

-1.94006 

.185872 

-.756482 

.176257 

-.427232 

.3125 

.209127 

-1.66725 

.3750 

.239824 

-2.66616 

.230487 

-2.44152 

.4375 

.239298 

-2.84153 

.5000 

.259715 

-3.50934 

.243199 

-3.06997 

.235401 

-2.91838 

.5625 

.219446 

-2.70160 

.6250 

.198364 

-2.30707 

.192882 

-2.20784 

.6875 

.157971 

-1.44573 

.7500 

.127328 

-.718524 

.120795 

-.468820 

.117826 

-.418792 

.8125 

.0764653 

.872537 

.8750 

.0397549 

2.43183 

.0388672 

2.42922 

.9375 

.0110272 

4.25203 

1 

0 

6.41809 

0 

6.40323 

0 

6.33992 
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TABLE  5  -  EXAMPLE  3:  THE  ADAPTIVE  MESH 


H 

64 

208 

634 

*1 

ujmgjj 

S32EH 

0 

0 

14.0294 

0 

20.5741 

0 

29.0727 

.0625 

.0366943 

9.02608 

.1250 

.0918232 

3.21720 

.0857160 

4.39572 

.1875 

.134002 

1.44502 

.2500 

.199010 

-1.94006 

.185894 

-.760749 

.176321 

-.433832 

.3125 

.209179 

-1.67517 

.3750 

.239801 

-2.67162 

.230509 

-2.45064 

.4375 

.239270 

-2. 8517 o 

.5000 

.259715 

-3.50934 

.243121 

-3.07583 

.235305 

-2.92978 

.5625 

.219247 

-2.71465 

.6250 

.198263 

-2.31303 

.192620 

-2.21778 

.7500 

.127328 

-.718524 

.120724 

-.475214 

.117593- 

-.429559 

.8125 

.0762835 

.858996 

.8750 

.039362 

2.42498 

.0387926 

2.41697 

.9375 

.0110087 

4.24102 

1 

0 

6.41809 

imam 

6.40290 

0 

6.33114 
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(4)  The  exact  value  of  u^x^.O)  (when  available) 

(5)  The  exaot  value  of  u^x^O)  (tdien  available) 

(6)  Number  of  elements  N  on  the  full  donain  (without  using  symmetry) 

Table  2  shows  the  results  for  Ex an pie  1.  In  this  case  the  exact  solu¬ 
tion  is  knoiei.  Good  results  have  been  obtained  especially  with  respect  to  the 
value  of  u1  (the  displacement).  The  accuracy  of  the  function  ug  is  also  good, 
but  not  quite  as  good  as  that  of  Because  the  solution  is  smooth,  the 
adaptive  mesh  is  uniform. 

Table  3  shows  the  results  for  Example  2.  The  solution  is  smooth  so  the 

adaptively  constructed  mesh  is  uniform.  The  exact  solution  is  not  available, 

but  the  results  improve  as  the  number  of  elements  increases.  As  before,  the 
function  u^  is  more  accurate  than  the  function  ug.  This  is  to  be  expected 
since  u2  =  Au^. 

Table  4  shows  the  results  for  Example  3  when  a  uniform  non-adaptive 
mesh  is  used.  The  exact  solution  has  a  strong  singularity  at  the  origin  and 
u2(x.|,0)  +  ®  as  *1  +  0,  The  exact  solution  is  not  available.  Once  more  the 
results  are  good  especially  with  respect  to  u^. 

Table  5  gives  the  results  for  Example  3  when  an  adaptive  mesh  is  used . 

Again  the  method  provides  good  results.  The  accuracy  for  an  adaptive  mesh  of 

634  elements  is  comparable  to  the  accuracy  provided  by  a  uniform  mesh  of  1024 
elements. 

The  results  displayed  in  these  tables  are  typical .  Other  test  computa¬ 
tions  were  performed  which  tend  to  confirm  the  general  behavior  shown  in  these 
tables . 


CONCLUSIONS 

Although  convergence  of  the  splitting  method  is  theoretically  guaranteed 
at  present  only  under  fairly  stringent  assumptions,  it  seems  to  work  very  well 
also  when  these  assunptions  are  not  satisfied.  The  method  is  a  mixed  one  of 
the  "saddle  point"  type.  Its  performance  is  very  good  with  respect  to  the 
displacement  u1  and  slightly  worse  with  respect  to  u2  or  u^.  It  is  apparent 
that  the  method  is  reliable,  particularly  if  it  is  controlled  by  an  error 
estimator.  Accordingly  it  is  recommended  for  use. 
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A  reliable  extraction  of  the  moments  and  shear  forces  from  the  computed 
results  would  appear  to  require  postprocessing.  This  problem  should  be  ad¬ 
dressed.  It  is  hoped  that  the  convergence  of  the  method  will  be  theoretically 
established  under  less  stringent  assumptions  that  those  imposed  at  present 
and  will  apply  to  the  examples  tested  in  this  report.  The  splitting  method 
is  preferable  to  the  reduction  to  the  Lame1  elasticity  equations. 

The  method  based  on  the  reduction  to  the  Lame'  elasticity  equations  is 
very  well  founded  theoretically,  but  practically  has  some  essential  disadvant¬ 
ages  which  makes  it  less  attractive.  The  main  disadvantages  of  this  method 
are: 

.  The  need  for  a  particular  solution. 

.  The  loss  of  accuracy  and  increased  error  due  to  the 
subtraction  of  the  particular  solution. 

.  The  necessity  of  assuning  that  the  domain  is  simply 

connected.  (This  restriction  simplifies  the  computational 
effort  considerably.) 

.  The  method  does  not  directly  provide  the  displacements. 

It  would  appear  that  the  splitting  method  is  the  preferred  method  for  the 


FEARS  user. 
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